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We study Gutzwiller-projected variational wavefunctions for charged, spinless holon excitations in 
chiral spin liquids. We find that these states describe anyons, with a statistical phase $s that is 
continuously adjustable between and it/2, depending on a variational parameter. The statistical 
flux attached to each holon is localized to within a lattice constant. By diagonalizing the effective 
Hamiltonian for charge motion in our variational basis we obtain a two-holon pair state. This 
two-holon state has large overlap with a two-hole state whose relative wavefunction has d-wave 
symmetry and breaks time-reversal and parity. 
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I. INTRODUCTION 



The notion of a quasiparticlc is a remarkably successful 
organizing principle for understanding many-body sys- 
tems. In strongly interacting Fermi systems such as "^He, 
for example, many-body effects simply renormalize the 
mass and interactions of the Landau quasiparticles; their 
charge_spin, and statistics the same as those of the bare 
atoms.EJ Other many-body systems, however, possess low- 
energy charged excitations that are dramatically differ- 
ent from Landau quasiparticles. In the fractional quan- 
tum Hall effect, for example, a large magnetic field and 
the Coulomb repulsion between electrons conspire to pro- 
duce Laughlin quasiparticles with fractional charge and 
statistics.B In trans-polyacetylene, a system with broken 
translational symmetry, charges bind to domain walls, 
formiiig soliton excitations with charge e but without 
spin.B In both of these cases, the novel quasiparticles have 
a "topological" character, because they cannot be cre- 
ated by local combinations of bare electron creation and 
destruction operators, in contrast with Landau quasipar- 
ticles. Nevertheless, these unconventional quasiparticles 
are well defined, weakly interacting objects which de- 
scribe the low-energy states of the many-body system. 

Similarly, unconventional quasiparticles have been sug- 
gested as the charge carriers in "resonating valen£&.h.ond" 
theories of high temperature superconductivityO'ErQ Mo- 
tivated by the fact that cuprate high temperature su- 
perconductivity arises upon doping antiferromagnets ex- 
hibiting large quantum fluctuations, these theories focus 
on spin- 1/2 systems with local antiferromagnetic corre- 
lations but no long range order - so-called "spin liquids." 

The phase diagram of the high temperature supercon- 
ductors motivates the study of spin liquids as a step- 
ping stone in a novel approach to the superconducting 
state. Referring to fig. |, the real materials reside in 



the temperature-doping plane; an additional axis, labeled 
frustration, is added for theoretical convenience. When 
the parent insulators are doped, antiferromagnetic long 
range order is destroyed almost simultaneously with the 
onset of superconducting order. (A "spin glass" regime 
due to disorder may intervene.) It is appealing to try to 
separate the effects of loss of magnetic order from the in- 
stability to superconductivity. In principle, this could be 
accomplished theoretically by first introducing fictitious 
frustrating interactions in the antiferromagnetic insula- 
tor, thereby destroying long range order without intro- 
ducing charge carriers. Then, to study the appearance 
of superconductivity, this "spin liquid" state would be 
doped, so that one could consider the motion of dilute 
charge carriers in a spin background resembling the fluc- 
tuating spins in the real materials' superconducting state. 

Finally, to return the theory to a more realistic model 
of the actual copper oxide planes, the fictitious frustrat- 
ing interactions would be turned off - they would be 
no longer needed, since magnetic frustration is implicitly 
generated by charge carrier motion.Q In this roundabout 
manner, one would arrive at a description of the real ma- 
terials' superconducting state. Assuming that no phase 
transition intervenes, the superconducting state achieved 
upon doping the spin liquid would have the same sym- 
metry and general features as the real superconducting 
state. This approach contrasts the more direct "BCS 
approach," which begins from an understanding of the 
normal metal at high temperature (which may be prob- 
lematic in the cuprates) and then considers residual at- 
tractive interactions between quasiparticles. If properly 
executed, both the spin liquid and the BCS approaches 
should arrive at the same superconducting state, but 
from different directions. 

Unfortunately for the spin liquid approach, it has 
proven difficult to find suitably frustrated spin mod- 
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FIG. 1. Schematic phase diagram for high temperature 
superconductivity. The real materials occupy the temper- 
ature-doping plane; an additional axis, labeled frustration, 
is included for theoretical convenience. In a BCS-like ap- 
proach, the superconducting state is entered from a well un- 
derstood normal state. In the "resonating-valence-bond" or 
"spin liquid" approach, a disordered quantum antiferromag- 
net is stabilized by artificially adding magnetic frustration 
to the undoped parent compound. This spin liquid is then 
doped to form a superconductor. The condensate is formed 
from bound holon pairs, the charge carrying quasiparticles 
in a doped spin-liquid. If no phase transition intervenes as 
the frustration is turned off in the doped state, both the spin 
liquid and BCS approaches must arrive at the same super- 
conducting state, but from complementary physical points of 



els which destroy antiferromagnetic order without in- 
troducing additional unwanted broken symmetry ground 
states. There is no short-ranged spin model known to 
have a spin liquid ground state. Nevertheless, spin liq- 
uids may be usefully thought of as metastable phases of 
frustrated antiferromagnets that mimic the spin corre- 
lations of cuprate superconductors, even if they are not 
the true ground state of a particular frustrated spin sys- 
tem. One can still study their excitation spectra, and 
hope to carry out the program outlined above. In this 
work |We consider simple variational "chiral spin liquid" 
stateaS which break time-reversal and parity symmetries 
but are translationally invariant. They have good vari- 
ational energies for simple frustrated spin Hamiltonians, 
and suggest natural variational states when doped. By 
studying the motion of a few charge carriers, we develop a 
scenario for understanding superconductors with broken 
time-reversal and parity symmetry. 



We consider the t-J-J' Hamiltonian H — H: 
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where the first term is the (frustrated) spin Hamiltonian 



(1) 



with first and second neighbor antiferromagnetic cou- 
plings J and J', and the second term is the nearest- 
neighbor hopping Hamiltonian 



cLc, 



(2) 



In these expressions c^^ creates a fermion at site i with 
spin projection a and Si = c\^<Tai3Cii3 is the spin operator 
at site i. Both Eqs. (|l|) and (|^) are allowed to act only 
within the subspace of states in which no lattice site is 
doubly occupied. 

In this work we address (a) the internal structure of 
charged "holon" excitations in a locally stable spin-liquid 
and (b) the dynamics and statistics of a dilute gas of these 
excitations. Sections |l| and HI contain a brief general dis- 
cussion of spin liquids and holon excitations and describe 
the specific wavefunctions we consider. Our main results 
are presented in sections IV and ^ A brief summary is 
provided in section VI. 



II. CHIRAL SPIN LIQUIDS 

Spin liquids are quantum antiferromagnets whose Neel 
order has been destroyed by large zero-point fluctuations 
of the interacting spins. These large quantum effects pre- 
clude the use of a spin- wave expansion about an ordered 
state. Instead, we appeal to the variational principle and 
take advantage of qualitative similarities between spin 
liquid states and certain simple Slater determinant states 
of spin-1/2 fermions. In particular, consider Slater deter- 
minant states Ix) which (a) are translationally invariant 
with an average density of one fermion per site, (b) have 
no average magnetic moment at each site, and (c) have 
spin and charge correlations which decay exponentially. 
All three of these properties are shared by spin liquids. 
The main difference between |x) and a spin liquid is sim- 
ply one of degree - a spin liquid has no number fluctua- 
tions at all, while |x) has local number fluctuations. 

This distinguishing property of \x) can be eliminated 
by applying the complete Gutzwiller projector 



(3) 



which annihilates all conflgurations containing doubly- 
occupied sites. What remains is a superposition of spin 
configurations whose amplitudes are inherited from the 
original "pre-projected" Slater determinant. If this pre- 
projected state has short-ranged spin and charge correL 
tions, then Gutzwiller projection is a "local" procedure 
The correlations in the projected state then closely follow 
the correlations of the pre-projected determinapt,|_t|p to 
simple renormalization factors of order unity.EjO^Ej 

It is surprisingly difficult to construct Slater determi- 
nants with the three properties listed above. To ob- 
tain a translationally invariant free-particle state with 
only short-range correlations, one should completely fill 
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FIG. 2. Schematic of a tight binding model which 
breaks time-reversal and parity symmetries. Circles repre- 
sent sites; lines represent hopping matrix elements between 
sites. A line with an arrow has phase 7r/4 in the sense given 
by the arrow, a dashed line has phase tt, and a solid line zero 
phase, (a) The flux state. Each elementary plaquette has vr 
flux piercing it. (b) The chiral state. Each elementary triangle 
has flux 7r/2 piercing it. Up to a gauge transformation, these 
tight-binding Hamiltonians are translationally invariant. 

a band of single-particle states while simultaneously en- 
suring that a gap exists between the filled band and the 
next highest energy band. For a translationally invariant 
system, the resulting charge density will be uniform. The 
difficulty arises in constructing a translationally invariant 
system with a one-body energy gap. The simplest way 
to accomplish this is to consider particles moving on a 
lattice in a uniform commensurate (fictitious) magnetic 
field which couples only to orbital motion.EHl3 For appro- 
priate field strength, the effective unit cell is doubled, 
and a gap is opened at half-filling. The hopping matrix 
elements will now be complex; as long as the net phase 
factor in the product of the elements around any closed 
loop is invariant under translations, the model ground 
state will also be translationally invariant, up to an over- 
all gauge transformation. 



FIG. 3. Single-particle bands of Eq. (|) in the two-site 
Square \/2 x %/2 magnetic unit cell shown by the dotted line 
in the inset, with no background magnetic field (dashed line) 
and in the presence of a translationally invariant chiral field 
(solid line) as illustrated in fig. |^ Note that a gap equal to 
4mo is opened at half-filling (filling the lower band). 

In particular, let |x) be the half-filled ground state of 
the square-lattice tight-binding HamiltonianEl 

Hx^^XtsclaCs-y (4) 

ija 

illustrated in fig. ^j. Here Xij is a complex link vari- 
able and creates a fermion at site i with spin pro- 
jection cr. We arrange for a net phase tt around ev- 
ery plaquette, and phase 7r/2 around every elementary 
triangle (encircled clockwise, say). More explicitly, the 
complex product XabXbcXcdXda around any elemen- 
tary plaquette ABCD has phase tt, while the product 
XabXbcXca around any elementary triangle ABC has 
phase 7r/2. We will refer to the phase around a closed 
loop a "flux," since the (Aharonov-Bohm) phase for mo- 
tion in a magnetic field is proportional to the flux en- 
closed by the path. One flux quantum hc/e corresponds 
to a phase 2tt. 

The Hamiltonian Eq. (^ explicitly breaks both time- 
reversal and parity symmetries, each of which inverts the 
flux through every closed loop. In this work we shall set 
Xij — for all but nearest and next-nearest neighboring 
sites i and j. We set \xij\ = 1 when sites i and j are 
nearest neighbors, and \xij\ = nio/'^ when sites i and 
j are next-nearest neighbors. This notation derives from 
the fact that for small values of mo and at half filling, the 
noninteracting fermions whose motion is governed by Eq. 
(^ have a gap to excitations equal to 4too ("mass," in 
field-theory jargon). Figure ^ illustrates the gap opening 
due to the addition of the diagonal hopping strength mo. 

The similarity of the Gutzwiller-projected states ^blx) 
to the actual ground states of a doped antiferromagnet 
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may be measured by their variational energy. Consider 
Ti with Jij = J for nearest neighbor sites i and j, and 
J'ik = J' for next-nearest neighbor sites i and k. Then 
for each value of J' / J there is a corresponding value of 
mo which minimizes the variational energy 



(7) 



E{J'/J;mo) 



{x\Pg-h,pUJ'/J)Pg\x) 



{x\Pg\x) 



(5) 



mo 



In this paper we study systems with values of toq sat- 
isfying < Too < 1, which implicitly selects the value 
of J' / J via Eq. (||) (see ref. |l^ for a more complete de- 
scription of the properties of E{J' / J; mo)). For small J' 
it is generally believedEJ that the ground state of Tispin 
has Neel order, unlike the variational states we consider 
here. For larger J' other more complicated types of bro- 
ken symmetry ground states may be stabilized. For our 
purposes, we assume that the spin liquid states Pg\x) 
are adequate models for the spin background of a lightly 
doped quantum antiferromagnet. 

The parameter toq plays a crucial role in our calcula- 
tion. It corresponds to the gap in the spectrum of at 
half-filling (see fig. ^) and determines the decay length 
for exponentially decaying spin and charge correlations in 
the pre-projected state |x). For small mo, ~ l/^o- 
(For general values of toq, ^_-£too) has a more compli- 
cated but nonincreasing formll3 when < mo < 1.) The 
spin-spin correlation length in the projected state Pfs]x) 
is related to by a numerical factor of order unity,E£l so 
Too also adjusts the range of spin correlations in the fully 
projected state. 

The introduction of a non-zero mo in Eq. (^) breaks 
both T (tjjpe-reversal) and P (two-dimensional parity) 
symmetryJj This can be seen, for example, by computing 
the spin triple product {Si-Sj x Sk) on any elementary tri- 
angle ijk in the projected state. This expectation value 
vanishes by symmetry when too vanishes, but is non-zero 
otherwise. 



III. SPINONS AND HOLONS 

Unlike the creation of quasiparticles in a Fermi liquid, 
the creation of separated spin and charge excitations in 
spin liquids involves non-local operations. There are two 
schemes for producing holons and spinous in Gutzwillcr 
projected states, which at first glance seem quite differ- 
ent. |— I 

Scheme 1 - AndersonllJ begins with a locally non- 
neutral state of the form 



Ci-<TcL|x), 



(6) 



a Slater determinant to which a particle with spin pro- 
jection cr' has been added at site j and a particle with 
spin projection —a has been removed from site i. In this 
state, {rij) — 3/2 and {rii) — 1/2, while (n^) = 1 for all 
other sites fc, as in the original Slater determinant state 
Ix)- The projected state 



has a uniform charge density everywhere. The act of pro- 
jecting Eq. (^) for far-separated sites i and j is therefore 
"non-local" in that density is transferred from site j to 
site i. In terms of spin, the state jicr, jcr') is an eigenstate 
of Sf and S^- with eigenvalue a and cr', respectively. The 
expectation value of spin at other sites is zero. These are 
"spinous." 

Making charged excitations without spin simply re- 
quires annihilating an electron at the center of each 
spinon. Anderson's state with "holons" at sites i and 
j is then 



(8) 



where A/"y is a normalization factor chosen so that 
{ijVf) — Sii'Sjj' + 6ij'6ji'. Note that the choice of a 
and a' in Eq. (g|) is arbitrary, apart from an overall phase 
factor. The holons defined by Eq. (||) are unconventional 
quasiparticles with novel statistics originating from the 
non-local Gutzwiller projection procedure. 

Scheme 2 - An alternate scheme for constructing 
spinons and holons begins from Slater determinants with, 
uniform charge density but non-uniform spin density.E£l 
In this approach the average charge density is one parti- 
cle per site even before projection, and the non-locality 
occurs in the making of the pre-projected Slater deter- 
minant. 

Spinons are constructed by considering tight binding 
models like Ti.^, but with additional tt fluxes localized at 
the spinon locations. The resulting Hamiltonian Ti^ is 
illustrated in fig. ^, and the energy spectrum with and 
without additional fluxes in fig. |5| Adding the localized 
fluxes leads to the appearance of midgap states as shown 
in fig. 1^ which are localized around the fluxes. Despite 
the absence of translational symmetry in Ti^, the— half- 
filled ground state \xd) has uniform charge densityjlj be- 
cause of particle-hole symmetry. (This contrasts with the 
non-uniform charge density in the pre-projected spinon 
state Ci-a-Cj^,\x) in scheme 1.) The spin density is non- 
uniform due to the singly-occupied midgap state, which 
is localized near the defective plaquette. As in scheme 1, 
spinons must be constructed in pairs, since the net flux 
through a periodic system must be an integer multiple 
of 27r. The spinon pair must be connected by a string 
of negated Xij 's to maintain a uniform field around each 
spinon (see fig. Essentially, the non-uniform charge 
density in the scheme 1 pre-projected spinon state (which 
leads to non-locality in the subsequent Gutzwiller projec- 
tion) has been replaced in scheme 2 by magnetic fiux in- 
sertion, which has long-range consequences for the tight 
binding matrix elements. 

Adding these localized fluxes is analogous to Laugh- 
lin's procedurea for creating quasiparticles in a fractional 
quantum Hall state, although here only half a, flux quan- 
tum is added. Since Gutzwiller projecting a locally 
charge-neutral state should not qualitatively change the 
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long-range correlations of the state, the extra tt flux at 
the defective plaquettes should survive Gutzwiller pro- 
jection, resulting in localized excitations with "semion" 
statistics halfway between fermion and boson. Charged 
holons are created in this scheme by removing a site at 
the center of the added localized flux. 

Despite scheme 2's intuitive appeal, we will not use it 
further, since we are interested in off-diagonal matrix ele- 
ments of the hopping Hamiltonian Eq. (0) between holon 
states at different positions. In scheme 2, this would re- 
quire working in a mixed basis of two Slater determinants 
that are ground states of different pre-projected Hamil- 
tonians. Remarkably, it can be shown for a special value 
of mo that schemes 1 and 2 are essentially equivalent 
(details are given in the appendix). In this work, all our 
calculations use scheme 1 (Eqs. ^, and||). 



IV. HOLON STATISTICS 



The charge motion in a spin liquid is governed by the 
hopping Hamiltonian Eq. The off-diagonal matrix 

elements of H in the holon variational basis contain 



FIG. 5. (a) Density of states resulting from the chiral flux 
Hamiltonian Tlx, El- (§) when mo = 1 (solid line). The 
dotted line spectrum was computed on a 512-site system with 
periodic boundary conditions and broadened with a width 
0.01. (b) Density of states resulting from the same 512-site 
system but with a pair of tt flux defects, TCd as illustrated 
in fig. ^ Fluxes have been placed as shown in the inset. 
Adding the fluxes leads to the formation of midgap states 
which split off from the upper and lower band, as shown. See 
the appendix for further discussion of these midgap states. 



FIG. 4. Schematic of defect Hamiltonian Hd used in 
scheme 2 (described in the text) for creating spinous and 
holons. (a) A localized flux is added to a plaquette, identified 
here by a small open circle. To do this, each link between 
sites crossed by the Dirac string (dashed line) originating at 
the added flux is multiplied by —1. At the other end of the 
string another localized flux is created (not shown), (b) A site 
is deleted from the system by cutting all links to it (removing 
all hopping matrix elements to surrounding sites). 



information regarding the kinetic energy and statistical 
phase of holons. In particular, these matrix elements are 
not simply given by the bare tij. The effective matrix 
element for one holon hopping from site i to site j in the 
presence of another holon fixed at site k is given by 

Tkij = -{jk\Hhop\ik), (9) 

where i and j are nearest-neighbor sites. 

We evaluate Eq. (||) using a method described 
previouslyEj Briefly, we consider a cluster of A^'r sites 
embedded in an infinite system. Correlations within the 
cluster are computed exactly, while correlations outside 
the cluster are accounted for only approximately. The 
thermodynamic limit is taken by considering clusters of 
various shapes and sizes centered on the sites i, j, and fc, 
and extrapolating to the infinite cluster limit. In prac- 
tice, this limit is reached when the cluster size exceeds the 
correlation length in the pre-projected Slater determi- 
nant Ix) . This method is substantially more efficient than 
exact projection of Slater determinants on finite toroidal 
systems (see below). 

The effective holon hopping matrix element T^ij con- 
tains information about holon statistics.cZI Computation 
of Tkij allows a direct check of the fractional statistics 
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invoked in the picture of Zou et aL.lli The hopping Hamil- 
tonian Eq. (||) restricted to the holon pair basis Eq. (H) 
is 



TIt 



ijk 



Tkij\kj){ki\. 



(10) 



In general, we can separate Tk ij into a magnitude and a 
phase, 



A(fc,l) • dl 



(11) 



where the phase is interpreted as the hne integral of an 
effective vector potential, A. This vector potential can 
be separatedli3 into two distinct parts: A(fc, 1) = Ao(l) + 
As(fc,l). The part of A which is independent of the 
position of site k is absorbed into Aq, and corresponds 
to a uniform background magnetic flux through which 
the holons move. The remaining part A5 depends on 
the relative holon position. 

Since the holon wavefunctions Ea_ (H) can be made 
symmetric under holon interchange, CJ the vector poten- 
tial As contains complete information about the holon 
statistics. Consider moving one holon around a part- 
ner which is held fixed at site k. The net phase factor 
for hopping one holon completely around another fixed 
holon is / As ■ dl = 2<i>5. When the moving holon path is 
made infinitely large, $5 is the statistical phase for holon 
interchange.tZl 

To implement this program, we first compute the back- 
ground vector potential Aq by calculating Tkij when k 
(the fixed holon location) is far from i and j (the initial 
and final moving holon locations). We find numerically 
that the background fiux is indistinguishable from the fic- 
titious uniform flux introduced in the pre-projected state. 
The vector potential Aq corresponds to a uniform flux 
of TT through each elementary plaquette, as explicitly in- 
cluded in Ti^ . To compute A5 we again calculated Tk ij 's 
on a path circling the fixed partner holon, and subtract 
the contribution from Aq. 

What is the "form factor" of the statistical fiux bound 
to each holon? One might expect that the bound flux 
would be distributed in a region of linear extent ^ about 
the holon. The statistics would only be well-defined for 
holons separated by this distance. This situation is sim- 
ilar to the case of the ^He atom, which is a composite 
of six fermions that behaves as a boson only on length 
scales longer than a Bohr radius, so that the internal 
atomic structure is essentially fixed. 

Surprisingly, we find that for the holon wavefunctions 
Eq.(||) the statistical flux around each holon is sharply 
peaked in its immediate vicinity, and that the radius of 
the bound flux tube is the lattice spacing, not the spin 
correlation length. Figure ^ shows the statistical field 
density As • d\ where the path C encloses one plaque- 
tte, plotted against the plaquette's geometrical distance 
from the origin. Most of the statistical flux is concen- 
trated in the core region immediately surrounding the 



FIG. 6. Statistical field density As ■ dl through the ele- 
mentary plaquettes surrounding a fixed holon, when mo — 1 
(filled circles), mo = 0.4 (open squares), and mo = 0.2 (open 
circles). Monte Carlo statistical error bars are approximately 
the same size as the circles. The inset illustrates the positions 
of the plaquettes whose field density has been plotted; each 
data point is accompanied by a letter corresponding to its 
location. 



holon (four plaquettes). Thus the phase shown in fig. 
represents the total statistical phase for far-separated 
holons. 

Figure ^ shows the statistical flux $5 threading a small 
2x2 square loop around a fixed holon as a function of 
nio- For small mg, the spin-spin correlation length tends 
to infinity, and the statistical flux through this small path 
vanishes. For mo of order unity, the spin-spin correlation 
length is essentially the lattice constant, and the statis- 
tical flux nears 7r/2. 

A slight negative contribution may be present outside 
the core (see fig. ^ especially the filled points repre- 
senting uiQ = 1). We are unable to determine whether 
this contribution is present for inter-holon separations 
exceeding ~ 5 lattice constants. Thus we cannot rule 
out the possibility that a small statistical fiux density at 
large distances from a holon integrates to give a signif- 
icant contribution to ^s- The existence of such a long- 
range tail to the flux distribution could make the holon 
statistics ill-defined, since the net phase would depend on 
the inter-holon separation. Neglecting this possible long- 
range part of the flux distribution, we conclude that the 
effective statistics of the excitations defined by Eq. (H) 
are continuously adjustable between values of zero (cor- 
responding to a boson) and 7r/2 (a "semion"). 

Using an analogy with the fractional quantum Hall 
effect, Laughlin and Zou have arguedE3 that the holon 
statistics in chiral spin liquids should be quantized as 
semions. Although the logic of this argument seems 
strong, our calculations (and the numerical calculations 
of ref. 18) find a statistical phase that is not simply 7r/2. 
Since here we have carefully considered extrapolation to 
the thermodynamic limit, this is not the source of the dis- 
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FIG. 7. Total phase accumulated by hopping one holon 
around another for the shaded loop shown in the inset, as a 
function of the mass parameter mo. This is the core contri- 
bution to twice the statistical phase, 2"l>s. The filled circles 
are computed using the cluster method, after extrapolation 
to the thermodynamic limit. The open circles represent the 
same quantity calculated on an 8 x 8 torus, which suffers from 
a discrete sampling effect and is inaccurate (see text). Monte 
Carlo statistical error bars are approximately the same size 
as the circles. The lines are guides to the eye. 



crepancy. There are two possible resolutions: (1) there 
is a weak, longer-range tail to the statistical flux distri- 
bution which we cannot resolve in our calculations, but 
which adds up to enforce <i>s = 7r/2, or (2) the variational 
states we have considered are not good holon states for 
small mo, but only for toq of order unity. It is only 
for such Too that the correspondence between scheme 
scheme 2, and the construction of Laughlin and Zouta 
has been firmly established. 

Our embedded cluster results agree with calculations 
performed on finite tori when properly extrapolated to 
the thermodynamic limit. While the two methods con- 
verge to the same result, the toroidal calculations con- 
verge much more slowly, and with important finitiL size 
effects. (This is also true for spin-spin correlations .113) To 
confirm this we calculated the core contribution to the 
statistical phase for L x L tori with L = 4, 6, 8, 10, 12, 
14, and 16. To confront the worst-case scenario for finite 
size effects we studied small toq in detail, where the cor- 
relation length ^ is longest. Fig. ^ shows the statistical 
phase from the nearest four plaquettes for toq = 0.01 and 
various toroidal systems. The embedded cluster method 
rapidly converges to zero statistical flux, as does the the 
sequence of tori with L = 4n -I- 2 = 6, 10, 14, ... The se- 
quence of tori with L = 4n = 4, 8, 12, 16, ... also extrap- 
olate to zero, but surprisingly slowly. This is likely due to 
the discrete Brillouin zone sampling at special symmetry 
points for these tori, a flaw which is not present in our 



FIG. 8. Extrapolation to the thermodynamic limit N —> oo 
of the core contribution to the statistical flux (shown in fig. 
^ when mo =0.01. (A^ is the number of sites in the projected 
Slater determinant.) Open squares reresent 4n x 4n toroidal 
systems, solid squares (4n+2) x (4n+2) toroidal systems. The 
4n X 4n toriodal systems suffer from an effect resulting from 
the discrete sum over the single particle Brillouin zone of H^, 
Eq. (^. Open circles indicate the result obtained with Ai'-site 
circular clusters embedded in an infinite system. The dashed 
line shows the result for a small embedded cluster (A^ = 13), 
at which point convergence has already been reached for the 
phase information in Tkij. The Monte Carlo statistical error 
bars are smaller than the markers. 

embedded cluster method. This accounts for the appar- 
ent discrepancy between the 8x8 torus results and the 
embedded cluster method in fig. |^. (The 4x4 torus is a 
particularly extreme case: the Slater determinant ground 
state wavefunction of is then independent of toq due 
to this sampling effect.) From this calculation we con- 
clude that the finite tori and embedded cluster methods 
agree when extrapolated to the infinite system limit. 

V. PAIR STATES 

We can now compute the pair wavefunction for two 
holons moving through a chiral spin liquid background 
by simply numerically diagonalizing Tix in the two-holon 
basis Eq. (^. We find that the lowest-energy two-holon 
bound state is a cprnplex d-wave state which breaks T 
and P symmetriesH3 

To diagonalize Tix, we note that this two-body Hamil- 
tonian is invariant under the same set of (magnetic) lat- 
tice translations as H^^. Bloch's theorem then states that 
each pair wavefunction can be labeled by a center of mass 
momentum K, and must have the form 

mi^^Y.e''''^'''^'''^^'<i>m^j), (12) 

where is the position of site i, and the relative wave- 
function (/)(ij) depends only on the lattice translation 
bringing site j into the same magnetic unit cell as site 
i, as well as the basis indices of the two sites. 
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TABLE I. Holon interactions. Hopping matrix element T 
for one holon fixed at site i and another moving from site j 
to site k when mo = 1. In the radial case sites i, j, and k are 
coUinear, in the azimuthal case the hopping direction (jk) is 
perpendicular to (ij). In each case djk > dij. 



dij 


\Tijk 

radial 


lA 

azimuthal 




1 


0.814 ±0.003 


0.827 ± 0.004 


1.20 ±0.08 


2 


0.806 ±0.002 


0.804 ± 0.004 




3 




0.811 ±0.002 


1.36 ±0.1 


5 


0.808 ± 0.002 


0.805 ± 0.003 


1.45 ±0.1 



By diagonalizing Tix we find the doubly charged 
ground state |<i>) in our variational basis. This calcu- 
lation of two holons in a chiral spin liquid resembles in 
spirit Cooper's original demonstration that a pairing in- 
stability exists for two electrons attracting above a Fermi 
sea. The effective holon hopping Hamiltonian Tix pro- 
vides a kinetic resonance energy of order a few percent 
t for nearby holons, because the magnitude |Tijfc| is en- 
hanced when sites i and j (or i and k) are close (see table 
^. Also, we have up until now neglected a static interac- 
tion Vij between holons at sites i and j which arises from 
the spin interaction, Eq. (|l|): 

= • ~ (^I^gS, • S™Pg|x». (13) 

(hn.) 

Roughly speaking, Vij provides a short-range static at- 
traction because neighboring holons can share a broken 
antiferromagnetic bond. In the undoped spin liquid, 
each unbroken bond contributes an energy J(S • S)nn = 
(—0.205 ± 0.005) J when toq = 1. To estimate the static 
holon potential we compute (zj|S/-S,„|ij) for the nearest- 
neighbor sites (Im) around a holon pair fixed at sites i 
and j using the embedded cluster method. Table ^ sum- 
marizes the results when mo = 1. A far-separated holon 
pair (5 lattice constants apart) breaks eight antiferro- 
magnetic bonds, costing 8 x |J(S • S)nn| ~ 1.6J, but 
the remaining nearby bonds are strengthened, so isolated 
holons end up costing slightly less energy to create. As 
the holons are brought closer together, they eventually 
share a broken bond, causing a net attractive potential 
of strength ~ 0.3J. Since both the kinetic resonance 
and the static potential are attractive, the variational 
ground state |<&)k is a bound state of two holons. Nu- 
merically, we find that the lowest energy bound state has 
center of mass momentum K = 0, and effective mass 
m* = l/{d^E/dK^) - 3m, where m ^ A\T\ is the ef- 
fective mass of an isolated holon with hopping matrix 
element T, evaluated for far-separated holons. 

What is the symmetry of the pair state |$)? To an- 
swer this question, consider the relative pair wavefunc- 
tion (piij). Since the relative phases of the holon pair 
states \ij) for different ij depend on the arbitrary choice 
of gauge used to construct the Slater determinant \x), the 
relative phase of ^(zj) are also gauge dependent. This 



FIG. 9. Overlap between (a singlet hole pair, 

Eq.(p^)) and \ij), (a holon pair, Eq.(^)) as a function of geo- 
metrical distance dij between sites i and j when (a) mo — 1, 
(b) mo = 0.4, (c) mo = 0.2. For a given pair, site j may 
lie on the same sublattice (open circles) or the opposite sub- 
lattice (filled circles) as site i. The correlation length is 
indicated in each frame. The dominant overlap occurs for 
nearest-neighbor pairs, and is close to unity over the entire 
range of mo. The overlaps are large for pairs closer than 
and are slightly larger for opposite sublattice pairs than for 
same sublattice pairs. The Monte Carlo statistical error bars 
are approximately the same size as the circles. 

complication can be eliminated by projecting the two- 
holon state |$) into the set of states comprising a singlet 
hole pair added to the undoped spin liquid: 

^ M^j (ctCjx - cxcji) PgIx)- (14) 

Here Mij is a normalization factor chosen so that 
((u) I (*'/)) — Sii'Sjj' + Sij'Sji'. These hole pair states 
are independent of the gauge choice in \x), since Pg\x) 
has precisely one particle per site, and is independent of 
gauge choice. The inner product 

m) = ((u)!*) (15) 
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FIG. 10. The relative pair wavelunction all (j) 

n are 

real numbers. The symmetry is a combination of d^2_y2 and 

is then a gauge invariant description of that part of the 
two-holon state which resembles a singlet hole pair. 

How similar are two holons and a singlet hole pair? 
The overlap between these two states is simply an expec- 
tation value in the state |x) 

{{ijm = 2M,X,,(x|Pg(1 - n,^)n,^4^c,i\x), (16) 

and can be computed using the embedded cluster me- 
thod: The magnitude of this overlap versus inter-holon 
distance dij is plotted in fig. ^ for three different val- 
ues of mo. A nearest- neighbor holon pair is essentially 
identical to a singlet hole pair, as shown by the near-unit 
overlap in fig. ^ at dij = 1 in all three cases. The overlap 
decays quickly as sites i and j are separated, presumably 
since far-separated holons disrupt the system differently 
than holes. Holes maintain a strong spacial correlation 
between spin and charge, while the holon has no spin. 
When holons are close together, however, the composite 
object has charge 2e and spin-0, permitting an overlap 
with a singlet hole pair with the same quantum numbers. 
The decay of the overlap with distance occurs on a length 
scale roughly equal to the correlation length as shown 
in fig. |. 

Using Eq. ( p^ ) and the two-holon state | $) obtained by 
exact diagonalization of Ht, we can compute the "hole 
pair wavelunction" (the inner product Eq. (|l5|)) for inter- 
holon separations up to 4 lattice constants. Figure [ll| 
shows the radial probability density as a function 

of the inter-holon distance dij. The amplitude is mostly 
concentrated near the origin, due to the kinetic resonance 
effect already mentioned and the short range of the holon 
pair/hole pair overlap. Further static local attraction due 
to the sharing of broken bonds only enhances this effect. 



FIG. 11. Probability density in the two-holon wavefunction 
$(y) versus inter-holon separation for holons governed by 
TCt with mo = 1 (filled circles), mo = 0.2 (open circles). In 
each case the wavefunction's angular symmetry is d-wave as 
illustrated in fig. 

without changing the angular symmetry of the state, as 
we verified numerically by considering additional short- 
range static attraction mimicking Vij . 

For all values of mo, the relative hole pair wavefunc- 
tion <&(ij) has d-wave symmetry (i.e., the state cha nges 
sign under 90 degree rotation), as illustrated in fig. |l^. 
Because of the explicit breaking of time-reversal and par- 
ity in the chiral spin liquid, this state is a complex linear 
combination of dx-2_y-2 and dxy states. (Although these 
states transform as different irreducible representations 
of the tetragonal group, they are mixed in the group of 
four-fold rotations that remains when reflections arc re- 
moved from the tetragonal group.E3) As shown in fig. |l^, 
the relative pair wavefunction is purely real on the x and 
y axes, and purely imaginary along the diagonals; '^{ij) 
is permitted by symmetry to be a general complex num- 
ber elsewhere, but we find that it remains either purely 
real or imaginary in non-symmetry directions as well. 

Since the two-holon state |$) has substantial overlap 
with a two-hole state, it is natural to consider a multi- 
holon state obtained by macroscopically occupying the 
hole pair wavefunction ^{ij), thus arriving at a BCS su- 
perconductor with a d-wave Cooper wavefunction. To 
show that this is the ground state of 7i or a similar spin 
Hamiltonian requires a much more elaborate calculation 
than we have given here. Nevertheless, if we assume that 
the "anyon condensate" of Laughlin et al. is obtained by 
a Bose condensation of holon pairs, then we may infer 
from our calculation that such a state will have conven- 
tional electron-pair off-diagonal long-range order: 

hm {cl^clc,,icj,^)^^*{ij)^{z'f). (17) 

The order parameter ^(ij) is nothing more than the in- 
ternal Cooper pair wavefunction, which for dilute, non- 
overlapping pairs should be given simply by the relative 
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hole-pair wavefunction ^(ij) we have calculated above. 
Thus a holon condensate is expected to possess all of the 
features of a BCS superconductor that has the uncon- 
ventional d-wave order parameter described above. Most 
strikingly, such a c?-wave state is expected to have a full 
energy gap, since there are no nodes in the Cooper pair 
wavefunction. 
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VI. CONCLUSIONS 



APPENDIX: EQUIVALENCE OF HOLON 
SCHEMES 



We have investigated the properties of charged excita- 
tions in a chiral spin liquid on a square lattice, using a 
variational basis which explicitly separates the spin and 
charge degrees of freedom. By directly computing the off- 
diagonal Hamiltonian matrix elements of the holon exci- 
tations, we found that the statistics of the excitations 
described by these variational states vary continuously 
between boson ($5 = 0) and semion ($s = 7r/2). The 
gauge field responsible for the statistical phase is local- 
ized in the area surrounding the holon, with a core size 
comparable to the lattice spacing. 

These results may be intuitively understood by view- 
ing holons and spinous aS|-tppological excitations in a 
chiral tight-binding modcl.E3 In the appendix we show 
that this picture (while computationally more difficult 
to implement) is an essentially equivalent representation 
of Andersons holons for spin liquids with short spin-spin 
correlation lengths. 

The two-holon ground state has total momentum 
zero, and its relative wavefunction transforms as a lin- 
ear combination of the d^2_y2 and dxy representations 
of the tetpagonal group. Such a linear combination is 
permittedE2l by the explicit breaking of T and P symme- 
try in the chiral spin liquid. Two factors contribute to 
the formation of a two-holon bound state: (a) a static 
attraction of order J between holons, and (b) a kinetic 
resonance due to enhanced holon hopping in the vicin- 
ity of another holon.L3 These effects are reminiscent of 
a "spin bag," in that the local environment around one 
holon lowers the energy of a second in its vicinity. 

The pair wavefunction of a bound holon pair has sig- 
nificant overlap with a bound singlet hole pair, espe- 
cially at short distances. This overlap suggests that a 
holon pair condensate is closely related to a hole pair 
(or electron) condensate, as found in the BCS theory 
of superconductivity. We conclude that a superconduc- 
tor arising from a pair condensate of bound holons is 
a BCS-like superconduptpr with a have a gap function 
with d-wave symmetry.Ej Such a superconductor would 
be fully gapped, since by mixing to id^y the gap func- 
tion avoids crossing zero. This superconductor can be 
described either as a d-wave BCS superconductor or as 
an anyonic superconductor. These two languages are ev- 
idently equivalent ways to describe the same state. 



In this appendix we show that the Anderson-holon and 
the topological charged excitations described previously 
in ref. |l^ are equivalent. This equivalence provides an 
intuitive picture of the Anderson-holon excitation and 
some insight into its local structure. The argument relies 
heavily on the equivalence between an Anderson-holon 
wavefunction and the KalmeyerpLaughlin wavefunction 
established by Laughlin and Zoulla so we start by briefly 
recalling the content of their proof. 

In ref. |l^ it is shown that the state obtained by 
Gutzwiller projecting the ground state of at half- 
filling and with toq = 2e^'^/'^ kJ1.912 is identical to the 
Kalmeyer-Laughlin (KL) state.Efl In first quantization, 
the KL state is expressed by the wavefunction 



N/2 N/2 



Zl, )^ Gizj) exp 
7=1 



\z,r 



, (Al) 



is a complex number defining 



where 

the position of an up spin labeled by 7 at square lattice 
coordinates (m,n), and G{z) is a gauge-dependent func- 
tion. G{z) — (_i)("+™+"™) in the gauge chosen by the 
authors of ref. |l8| . 

The proof consists of two steps. In the first it is shown 
that given the first-quantized version of the up spin part 
of the half-filled ground state of - a function 0({r}) 
where the set of coordinates {r} = ri, . . . , describes 
the position of N/2 up electrons on a lattice of N sites 
- one can write down the wavefunction corresponding 
to the Gutzwiller projected state ^' as a function of the 
same set of coordinates: 



^{{r}) = ^'{{r})e 



»/({r}) 



(A2) 



where the function /({r}) depends on the c hoice of the 
basis used to represent 5'. (See Eq. (A12) for a pre- 
cise definition.) In-Uae second step it is shown that with 
a particular choicecj of niQ, the Slater determinant ob- 
tained by filling the lower band of coincides with the 
determinant obtained by filling the energy levels corre- 
sponding to an electron in free space moving in a uniform 
magnetic field with Z = (Hc/eB)^^^ = provided the 

single-particle wavefunctions defined on a continuum are 
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restricted to points on a square lattice with lattice con- 
stant equal to unity. This representation of 4>, together 
with Eq. ( |A^ is then used to show that ^' = ipK l ■ As a 
simple corollary, Zou and Laughlin show that the wave- 
function corresponding to an Anderson-holon pair at sites 
A and B is represented by 



\za,zb) = Y\_{z-f-ZA)(zj - zb) W_{ztj. - z„y 



lJ.<v 



z,y ) exp 



(A3) 



Our approach is very similar. First we show that the 
first step follows generally when a particle-hole symme- 
try operator exists, thus it also applies to the class of 
Hamiltonians Hd described in ref. |l6|. Second, we argue 
that a representation of the filled lower band of Hd is 
obtained by filling the first Landau level corresponding 
to an electron in a uniform field with two flux tubes of 
appropriate strength at position ta and tb, again having 
restricted the wavefunction to points on a square lattice. 
Although we will not prove this point rigorously it is a 
natural assumption in view of the results of ref. 

The Hamiltonian Hd describing the topological 
charged excitation of section ||'s scheme 2 is obtained 
from H^ by adding fiux defects in the form of an addi- 
tional half-fiux quantum tt through a plaquette adjacent 
to site A, the same through a plaquette adjacent to site 
B, and a string of tt phase additions to the phases of the 
Xij's between sites A and^S to restore the proper fiux in 
the remaining plaquettes.E3 (The actual position of this 
string is irrelevant, as it can be moved by a suitable gauge 
transformation.) Note that the added fictitious fiux cou- 
ples only to the particles' orbital motion. 

It is convenient to start with H^. (and Hd) expressed 
in a gauge such that the Xij are real for nearest-neighbor 
and imaginary for next-nearest neighbor bonds. This 
choice of phases for the link variables allows a simple def- 
inition of the particle-hole conjugation operator T, which 
we define to be the linear operator satisfying 



TlvaCn 



(A4) 



where a is an arbitrary phase, Q = (tt, tt), jvaCo-) is 
the vacuum state in the a spin subspace, and \Fa) = 
J^^j.-j cJ|vaCo-) is the state with one a electron at each 
of the N sites in the system. To avoid ordering am- 
biguities this product is performed in an (arbitrary) 
order specified by an ordering function 0(r) number- 
ing the sites of the lattice from to — 1 such that 
0(ri) <0(r2) <---<0{vm). 

With this definition of T one can show that T is uni- 
tary and = e*"*", where 7 depends only on N and 
a. The usefulness of T is apparent when we realize 
that T is a symmetry operation for both Ti.^ and T-Ld'- 
[T, H^ — [T, Hd] = 0. Notice that both charge density 



and current density are odd under particle-hole symme- 
try, that is T^(l — nr)T = nr — 1, JaT ~ — Ja, where 

•/a = »E.(e'*^4.Cr+a,.-H.C.). 

We also define the antilinear operator Q : 



Qlvaccr) 

5 ^raQ 



\F„ 



(A5) 



Q is antiunitary, anticommutes with H^ and Hd, and 
commutes with T. A simple consequence of the commu- 
tation relations of Q, T, and the Hamiltonians Hd or H^ 
is that single-particle states always come in pairs, \E^) 
and QT\Ec), corresponding to energies E and —E re- 
spectively. Since in the absence of fiux tubes and at half 
filling the ground state of H^ is fully gapped, syjametry 
suggests (and calculations have explicitly shownllj) that 
with a single fiux-tube and an odd number of sites the 
energy spectrum of Hd will have a midgap state, whose 
amplitude is localized near the flux tube. With an even 
number of sites and two fiux-tubes at ta and tb the de- 
generacy of the two midgap states in the energy spectrum 
of Hd is lifted by A_Eioc, which is exponentially small in 
Ir^-rsl. (See fig. |.) 

Our task now is to evaluate the Gutzwiller-projected 
two-holon wavefunction of ref. This is obtained by 
first considering the ground state at half-filling \xd) of 
Hd when fiuxes have been added at sites A and B, 
and these same sites have been removed from the lat- 
tice. Notice that the removal of the two sites eliminates 
the two midgap states and that the ground state is non- 
degenerate and hence invariant under T. 

After Gutzwiller projection every configuration has a 
single electron at each site and an equal number of up and 
down electrons. A convenient basis in which to represent 
PclXd) is 



IW) 



(A6) 



where — c]!j^Cr| is the usual spin lowering operator, 
{r} denotes the set of positions of the up spin electrons, 
and {f } denotes the set complementary to {r} specifying 
the positions of the down spin electrons. The sets {r} 
and {f} each contain M = N/2 — 1 members. 

The Gutzwiller-projected wavefunction in this basis is 
given by 



vl/({r}) ^ iMlPclxd) = ({rllx.), 



(A7) 



since none of the states |{r}) contain site double- 
occupancy. Using the anticommutation properties of 
fermion operators, we find 



({r}|Xd) = exp 
where 



M 



M{r})M{f}), (A8) 



\Xdf/ 



(A9) 
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and we have made the decomposition \xd) = Ixdt)'^ IXdj.) 
and dropped the exphcit spin indices on Cr and cj. We 
can now use the properties of T to obtain 

M{i}) = (vact|cf„ ■■■Cf, Ixdt) 

= (vact|TTtcf„T...Ttcf,TTt|xdT) 

M 



exp 



exp 



1=1 

M 



•••CfJXdT) 



(AlO) 



i^Q-f, + 7rO(f,) 



(vaCtlCr^, •••CrilXdr; 



which, together with Eq. (A8), gives 



*({r}) = exp 



M 



z^Q 



(All) 



In Eq. jAll ) we have used the freedom in the choice of a 
in Eq. ( |A4D to chminate an additional phase independent 
of {r}. 

Eq. ( IaTiI ) gives the wavefunction 5* in the particular 
gauge in which the link variables Xij either real or 
purely imaginary. Under any local gauge transformation 
Crcr exp[— iA(r)]cr(T wc finally derive, with the help of 
Eq. (P), 



v|/({r}) = exp 



M 



2A(r,) 



0d({r};A), 



(A12) 



where we have added explicit dependence on A to 0^ in 
order to remind the reader that the up spin part o f the 
ground state \xd^) is gauge dependent. Eq. (A12) now 
corresponds directly to Eq. (A2), as desired. 

Next, we want to find an alternate representation of 
(j^d based on a continuum picture of electrons moving in 
a uniform magnetic field with two flux tubes at posit ion s 
A and B. Notice that although by its definition Eq. (A7) 
^' is gauge invariant, (pd is not, and we must be careful 
to choose a gauge in the continuum which corresponds 
t o the particular gauge used to derive Eq. (All) and Eq. 
(|A12|). In the following we work in the symmetric gauge 



although the gauge used to derive Eq. (All) and used in 
ref. m to express \za, zb) is the Landau gauge. (Strictly 
speaking, this gauge differs from the Landau gauge by a 
"singular" transformation, which will be discussed later.) 
We start with the continuum Hamiltonian 



H 



^' (-zV + A)' 



2toP 



(A13) 



where A = (r x z)/2 + Ar^ + Ar^, the unit of length 
I = {hc/eBy^'^ , and B has been chosen so that in the 
absence of flux tubes the degeneracy of the Landau levels 
is equal to N/2; this implies B = (j)Q/2 where <j)Q = hc/e 
is the flux quantum. The vector potential corresponding 
to a flux tube of strength a(f)Q is then given by 



(A14) 



where 9a is the azimuthal angle measured with respect 
to r^, and positive values of a correspond to a locahzed 
flux which strengthens the background magnetic field. 

It is convenient to gauge away Ar^ and Ar^ by looking 
for a solution in the form ijj — ijj exp[— i J(Ar^ +ArB)-(il]. 
Then ?/; obeys the Schrodingcr equation for a particle in 
an uniform field but must be multivalued in order to make 
il) single- valued: when encircling ya and yb the phase of 
ip must change by ua and as respectively, in order to 
compensate for t he ph ase change due to the exponen- 
tial term in Eq. ( A15 ). But the most general solution 
corresponding to the lowest Landau level orbitals in the 
symmetric gauge can be written as Tp^z) = g(z)e^^^^, 
where z — x + iy and g{z) is an arbitrary function of 
z. The solution of interest is then given by g{z) = 
{z — za)°''^{z — zb)"'^ g{z), with g{z) analytic everA^where. 

Although the solution so defined is valid for anynJ value 
of a A and a b , our purposes lead us to consider the spe- 
cial case a A = —cub = a > 0. In this case, as a is 
increased from zero, a lone state splits off from each of 
the n = and n = 1 Landau level manifolds: the first 
lone state's energy increases, and the second lone state's 
energy decreases. Note that the degenerate Landau level 
manifolds can be associated with the lower (n = 0) and 
upper {n — \) band of the lattice Hamiltonian Hd- The 
two lone states which split off are similarly associated 
with the two localized states of Hd in the vicinity of the 
flux defects. Filling the n = Landau level corresponds 
to computing the Slater determinant obtained by filling 
the levels 



xpkiz) = (z - zaT{z - zbY 



(A15) 



where fc 0, 1, • • •. The resulting Vandermonde determi- 
nant can be easily calculated giving the result 

</>(zi, • • ■ ,Zm) = ]J(2:-y-2:A)"(2;7 - ^s)^"" II (^/^ ~ 



l[Giz,) 



exp 



|2 " 



(A16) 



To complete the argument we now employ the gauge 
freedom associated with A in Eq. ( A12). First notice that 



the use of -ijj in place of tp to obtain Eq. ( A16) is justified 
because the gauge we were effectively using in writing a 
lattice Hamiltonian with a string of defects connecting 
sites A and B is precisely the singular gauge (singular 
only in the continuum limit, of course) we employed to 
go from 'ijj to ijj. Alternatively we could have used ip to 
define 4> in Eq. ( A16) and performed a gauge transforma- 
tion on the lattice Hamiltonian to smooth out the string 
connecting A and B. When this gauge transformat ion i s 
properly taken into account as prescribed by Eq. (A12) 
the result for ^ is identical. Finally we have to convert 
the Landau gauge, corresponding to Eq. (All), into the 
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symmetric gauge used to express 4> in Eq. (A16). This 
is easily accomplished through the gauge transformation 
defined by A(r) = TTxy/2 which, when restri cted to the 
lattice r = (m, n), makes the prefactor in Eq. ( A12 ) equal 
to 



exp 



M 



i^Q-r, + 2A(r, 



(A17) 



This is identical to the factor G{zj) appearing in Eq. 



By using Eq. (A12) and Eq. (A16) with the continuum 



variables restricted to lattice points and a — 1/2, 
the equivalence between the two-holon state defined as 
in Section II, scheme 2, a nd the Anderson holon state, 
scheme 1, defined by Eq. (A3) is now explicitly shown. 
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